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Abstract 

In this paper we investigate decoherence time of superconducting Josephson charge qubit (JCQ). Two kinds of methods, 
iterative tensor multiplication (ITM) method derived from the qusiadiabatic propagator path integral (QUAPI) and Bloch 
equations method are used. Using the non-Markovian ITM method we correct the decoherence time predicted by Bloch 
equations method. By comparing the exact theoretical result with the experimental data we suggest that the Ohmic noise 
plays the central role to the decoherence of the JCQ. 
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INTRODUCTION 



Solid state qubits are considered promising candidates 
for making processors of quantum computers because 
they can be scaled up to a large numbers. The qubits 
based on Josephson junction are these kinds of qubits. 
But how about their other qualities, in particular, how 
about their coherence? Many efforts not only theoretical 
0, but also experimental E [H El have been 
contributed to search decoherence time as well as deco- 
herence mechanism of the qubit systems. The theoretical 
researches are in general based on the spin-boson model 
By now, it is suggested that there are 1// noise and 
Ohmic noise in the environment of the Josephson qubits. 
It is considered that the former is derived from the back- 
ground charge fluctuations and the latter from the elec- 
tromagnetic noise due to voltage fluctuations. But what 
is the most primary mechanism of the decoherence, or 
which is the mainly noise source in the environment of 
the qubits? By using the Bloch equations [f| one can 
derive the relaxation time and dephasing time of the 
qubits. However, in the derivation of Bloch equations 
an approximation scheme in general the Markov approx- 
imation should be used. It has been pointed out recently, 
the Markov approximation is not a suitable approxima- 
tion scheme in the investigations of the quantum system 
for quantum computing purposes because it is not usu- 
ally valid at low temperatures and for short cycle time of 
quantum computation |jj . A similar viewpoint on a qubit 
of double quantum dots is also pointed out, see as @- 
Privman and his co-worker investigated the decoher- 
ence of qubits with short-time approximation rather than 
the Markov one, and many interesting and novel results 
are obtained. However, it is not enough to only investi- 
gate the decoherence in a short time because the coher- 
ence in a longer time for some qubits (for example qubits 
for quantum registers) is also interesting. Fortunately, 
it is found out that by using short-time propagators one 
can construct a path integral |lfj| in a longer time. The 
well established iterative tensor multiplication (ITM) al- 
gorithm derived from the qusiadiabatic propagator path 
integral (QUAPI) H3 can be used to solve the evolutions 



of low-dimension open quantum systems in a moderate 
long time 12]. In this method the temporal non-local 
interactions is involved and it is non-Markovian. Thus, 
we expect that the non-Markovian method can be used 
to investigate the decoherence of Josephson qubits and 
then to detect the mainly mechanism of the decoherence 
in the qubi ts. There are many kinds of Josephson qubit 
models [H Q E Q [H E [|! HJ , but in this paper we 
only investigate the Josephson charge qubit (JCQ). We 
shall obtain an accurate decoherence time of the JCQ by 
using the ITM algorithm without the Markovian approx- 
imation. Base on the accurate decoherence time we shall 
suggest a mainly mechanism of the decoherence in the 
JCQ. 



II. MODELS 

The elementary unit of the quantum computer is a 
qubit which is in fact a two-level quantum system [2l| . 
There are many physical realizations for the system. 
However, any single qubit can be represented as a spin- 
1/2 particle, and its Hamiltonian can be written as 
H (t) — —T;B(t) ■ a. Here o~ x . V:Z are Pauli matrixes in 
a space of states |t) = (0 1) and ||) = (l 0) which 

form basis states of the Hamiltonian. The quantity B (t) 
has different physical meaning according to the difference 
of the physical realizations of the qubits. For example, if 
the qubit is realized by a spin of some particle the B (t) 
will be an effective magnetic field. But to the JCQ, the 
components of the "magnetic field" are set as B x = Ej, 
where Ej is the Josephson energy of the Josephson junc- 
tion, B y = 0, and B z = AE C (1 - 2» 9 ) @. The Hamilto- 
nian of a general qubit can be represented as 



H s — — -B z a z — -B x a x . 



(1) 



If one modulates the gate voltage and makes n g = 1/2, 
the JCQ system may be reduced to H s = — ^B x a x . In 
this paper, we take Ej — 51.8 fieV, and Ec = 122 ixeV 
according to Ref. [2^ ■ If one considers the interaction of 
the qubit and its environment, and takes the environment 



1 



as a bath, the whole Hamiltonian of the system-bath will 
be 



where 



H = Hq + H e 



(2) 



Xk 



(3) 

In general, Hq should plus counterterms — J2 k \/^ rn k UJ k 
which can ensure that some important features of the 
qubit do not depend on the coupling strength. In our 
problem, the counterterms only contribute a global phase 
so we can ignore it. This is the well-known spin-boson 
model, a appropriately reduced open qubit model. In the 
following, we firstly analyze this model. On the one hand, 
for the bath, only the linearly coupling term is chosen 
in the coordinates Xk, representing the lowest nontrivial 
term in the Taylor series expansion of the potential. It is 
accurate enough in the weak coupling case. On the other 
hand, for the qubit, only a z coupling term is included. 
The terms of a XyV coupling with the bath have not been 
included in the Hamiltonian. The reason is that a XlV 
have only off-diagonal matrix elements in the a z repre- 
sentation, i.e., they only change |f ) to j) and vice versa. 

In order to obtain the reduced density matrix of the 
qubit in qubit-bath system, one should know the coupling 
coefficients Xk in Eq.|0). However, we do not need to 
know their details because all characteristics of the bath 
pertaining to the dynamics of the observable system are 
captured in the spectral density function of noise 



A 



m k LOl 



(4) 



In the case of truly macroscopic environment the spectral 
density is for all practical purposes a continuous function 
of frequency. In the following we shall see that the spec- 
tral density function J (u>) instead of Xk is directly used 
in obtaining the elements of the reduced density matrix. 
It is related to the power spectrum of the noise as 



S x M = Jx (uj)hcoth (ujf3h/2) . 



(5) 



Here, (3 = 1/ksT, where ks is the Boltzmann constant 
and T is the temperature. Throughout the paper we take 
T = 30 mK according to Ref. |22|. Due to very wide of 
the real and imaginary parts of the response function (see 
Ea. (|14|l in the following) in the time range for the 1// 
bath, we in fact cannot investigate the evolutions of the 
reduced density matrix of qubit in the 1/ / bath with the 
ITM. Thus, in this paper, we focus on the case that the 
environment is the Ohmic bath. The spectral density of 
the Ohmic bath can be expressed as [6j 



(6) 



Jx (w) = 2nhauj exp 



where u>c is the cut-off frequency of the bath modes. The 
parameter a is dimensionless strength of the dissipation 
which is determined by concrete qubit-bath system. For 
the JCQ model, Makhlin et al. proposed a numerical 
simulation value a « 10~ 6 . In this paper we suppose 
a = 5 x 10~ 6 . If a > 5 x 1CP 6 the decoherence time will 
be shorter than the results obtained in this paper. On 
the other hand, if a < 5 x 10~ 6 the decoherence time will 
be longer than the results in this paper. 



III. QUAPI AND ITM 

In the following, we firstly review the QUAPI method 
and then the ITM scheme. Suppose the initial state 
of the qubit-bath system has the form 



R(0) = p(0)® Pbath (0) 



(7) 



where p (0) and Pb at h (0) are initial states of the qubit 
and bath. The evolution of its reduced density operator 
of the open qubit 

p (s", s'; t) = Tr fcQt , ( S "| e~ im l h p (0) ®p hath (0) e lHt ' h \s') 

(8) 

is given by 

p(s\s'-t) 
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JHgAt/H 
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where the influence functional is 



') s jV-l! S i At) j 

(9) 



J ( s ' s l ' ' ' '' ^N— U s ' s ' s l ' ' ' '' ^N— 1' s ' 



TV, 



bath 



S N — 1' s ' s ' S l ' ' ' '' S T 

iH en v(s")At/2H -iHe nv (s+_ 1 )At/2h 



x . . . e -«r.«« ( s o ) At / 2h Pbath (0) e ^(s-)At/2H 

X . . . e iHe n -v{s-_ 1 )AtJ2h iH tnv {sl)At/2H 



(10) 



The discrete path integral representation of the qubit 
density matrix contains temporal nonlocal terms 



i {4,4 



Sjy_ 1 , s'; At) which 



denotes the process being non-Markovian. With the 
quasiadiabatic discretization of the path integral, the 
influence functional, Eg. (|1U|I takes the form 



N 



I = exp\--Y, E (4 ~ s k) (Vkk'4> -V*kk'3 k >) 
^ k=0 k'=o 

s" and s 



(11) 

N — s'. The coefficients i] kk , 
can be obtained by substituting the discrete path into 



where 
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the Feynman- Vernon expression. Their expressions have 
been shown in Thus, the influence functional can 

be expressed with a product of terms corresponding to 
different Afc as 



where 



JV N—l N-Ak 

i = uio(4) n n iA k {s±,s± 

k=0 k=0 k=0 



k+AkJ 



N-Ak 

n 

k=0 



(12) 



Here, Ak = k — k', where k' and k are points of discrete 
path integral expressions, see Ref . . and 



lAk (sf, S f +Ak ) 



exp 



He 



r lii s i i 



1 

h 



X (Vi+Ak,i s i ~ Vi+Ak,i S i )} , Afc ^ 1. 

(13) 

The length of the memory of the time can be estimated 
by the following bath response function 



lit) 



1 



coth 



cos Lot — i sin wt 



(14) 

It is shown that when the real and imaginary parts be- 
have as the delta function 5 (t) and its derivative 8' it) , 
the dynamics of the reduced density matrix is Markovian. 
However, if the real and imaginary parts are broader than 
the delta function the dynamics is non-Markovian. The 
broader of the Re[7 (t)] and Im[7 (t)] are, the longer of 
the memory time will be. The broader of the Re[7 (t)] 
and Im[7 (t)] are, the more serious the Markov approxi- 
mation will distort the practical dynamics. In Fig.l we 
plot Re[7 (t)] and Im[7 (t)] of the Ohmic bath. 



Fig.l 



From Fig. 1 we see that the memory time is about 
Tmem = 1 x 10 -11 s for the Ohmic bath. Due to nonlocal- 
ity of the influence functional, it is impossible to calculate 
the reduced density matrix by Eq. J^J in matrix multipli- 
cation scheme. However, the short range nonlocality of 
the influence functional Eci. (|10f) implies that the affects 
of the nonlocality should drop off rapidly as the "interac- 
tion distance" increases. In the Makri's ITM scheme the 
interaction can be taken into account at each iteration 
step. The reduced density matrix at time t = N At (N 
even) is given as 

p (4, NAt) = A« ( S ±; NAt) I (s%) , (15) 



AM (* + l)At) = J ds±T™ ( S ± a ± +1 ) (16) 

x^W (sf;kAt) . 



Here, 



k-\-A.k —1 
n—k 

xI 2 (st, -■ / Afc max (4 > S n+Afc max ) . ( 17 ) 



and 



A (AW) ( s ± 



'1 1 



Afe„ 



,-i,o) = (4\p. (°)l s o) 



where 



(18) 



fe+i 



) = <- 



x /s 



+ +1 \exp(-iH At/h) |s+) 
k] eMiH At/h)\s- +1 ). (19) 



In the ITM scheme a short-time approximation instead of 
the Markov approximation is used. The approximation 
makes a error of short-time propagator in order (At) 3 
which is small enough as we set the time step At very 
small. It is shown that when the time step At is not 
larger than the characteristic time of the qubit system 
which can be calculated with H/Ej the calculation is ac- 
curate enough 0- In particular, the scheme do not dis- 
card the memory of the temporal evolution, which may 
appropriate to solve the decoherence of qubits. 



IV. DECOHERENCE OF JOSEPHSON QUBITS 

To measure effects of decoherence one can use the en- 
tropy, the first entropy, and many other measures, such 
as maximal deviation norm etc., see 0- However, es- 
sentially, the decoherence of a open quantum system is 
reflected through the decays of the off-diagonal coherent 
terms of its reduced density matrix. The decoherence is 
in general produced due to the interaction of the quantum 
system with other systems with a large number of degrees 
of freedom, for example the devices of measurement or 
environment. In this paper, we investigate the decoher- 
ence time of the JCQ via directly describing the evolution 
of the off-diagonal coherent terms, instead of using any 
measure of decoherence. In our following investigations, 
we suppose the cut-off frequency of the bath modes is 
wc = 5 (ps) . We set the initial state of the qubit 
p(0) = 5 (|0) + |1)) ((0| + (1|) which is a pure state and 
it has the maximum coherent terms, and the initial state 



of the environment p bath (0) = JIa 



-l3M k 



/Tr fc (e-^). 



Decoherence time obtained from ITM scheme: At first, 
we use the ITM scheme investigating the decoherence 
time of the Josephson qubits. The evolutions of the co- 
herent elements of the reduced density matrix of the JCQ 
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in the Ohmic bath is plotted in Fig. 2. Here, we simply 
choose Afc max = 1 and At = 1.27 x KT 11 s in the ITM 
scheme. The choice of the time step is feasible as we con- 
sider that it should be not shorter than the memory time 
of the bath, because the latter is about T mem = 1 x 10~ n 
s for the Ohmic bath, see Fig. 1. It is also appropriate as 
we consider that the time step should be not longer than 
the characteristic time of the qubit, where the latter is 
about r = 1.3 x 1CT 11 s. 



Fig.2 



It is shown that when we choose the parameter of the 
dimensionless strength of the dissipation a — 5 x 10~ 6 , 
the decoherence time of the JCQ is about T2 = 1.05299 
/is. 

Decoherence time calculated on Block equations: It is 
well known that the decoherence time can be derived 
on Bloch equations. In this method, the relaxation and 
dephasing times r% , T2 can be calculated as [23] 

rr 1 = 2r^ = J (w„) coth (/3^ /2) , (20) 

where loo = B x /H is the natural frequency of the Joseph- 
son qubit. From Eq. (|20|) and using the same parameters 
of the qubit and the bath as above we can obtain that 
the decoherence time is T2 = 1.61966 /is. It is shown that 
the time obtained from Ea. H20|) is longer than that from 
the ITM scheme. We suggest that the difference is de- 
rived from the following two reasons. The first is that the 
Bloch equations are in general derived from the Markov 
approximation which discards the memory of the bath in 
the derivation of the dynamical evolution. The second is 
that the Ea. H20|) is obtained from the second order ap- 
proximation of perturbation series. The decoherence of 
the qubit described with this method is only the "res- 
onant decoherence" [24|. It is not equals to the actual 
decoherence accurately except for the "nonresonant de- 
coherence" very small. 

Compared with the experimental results: In our calcu- 
lations, we use the parameters similar to Rcf. 22J, so we 
can compare our results to the experimental decoherence 
time. In the decoherence time of a single-Cooper pair 
box, namely, the JCQ is estimated. The main decoher- 



ence source is thought to be spontaneous photon emis- 
sion to the electromagnetic environment (which can just 
be described by the Ohmic bath). In [22] the authors 
pointed out that the experimental decoherence time of 
the JCQ could exceed 1 /xs. It is shown that by use of the 
Ohmic decoherence mechanism we can obtain a theoret- 
ical decoherence time of the JCQ not only by the ITM 
scheme but also through the Bloch equations method. 
Both of the theoretical results are agreement with the 
experimental one very well! 
V. CONCLUSIONS 

In this paper we investigated the decoherence time of 
the JCQ in the Ohmic bath with the ITM scheme based 
on the QUAPI and based on the Bloch equations. The 
results derived from the two kinds of methods are com- 
pared with each other. It is shown that the decoherence 
time obtained from the Bloch equations method is longer 
than that from the ITM scheme. We suggest that the dif- 
ference is resulting from the different choices of the ap- 
proximation scheme because the Markov approximation 
used in the Bloch equations method discards the mem- 
ory of the bath. It is also because the Bloch equations 
method discards the higher order decoherence, namely, 
only the "resonant decoherence" [24| is left over. So the 
decoherence time obtained from this method is not equals 
to the actual decoherence time accurately. What is more 
important to us is that the experimental decoherence 
time of the JCQ due to spontaneous photon emission 
is well agreement with the ITM decoherence time of the 
JCQ because of the electromagnetic fluctuations. Both 
of the spontaneous photon emission and the electromag- 
netic fluctuations have the same decoherence mechanism 
and can be modeled by the Ohmic bath. These can lead 
to a conclusion that the Ohmic bath decoherence is a 
central mechanism in the JCQ and the decoherence time 
of the JCQ is about 1 fis when the temperature is about 
30 mK and the Josephson energy is about 51.8 /xev. The 
decoherence time is decided by the decoherence mecha- 
nism and affected by the experimental temperature. If 
the experimental temperature increase the dimensionless 
strength of the dissipation a will also increase and the 
decoherence time will be shorter, and vice versa. 
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VI. CAPTIONS OF THE FIGURES 

Fig.l: Real (line) and imaginary (short line) part of 
the response function of the Ohmic bath. Here, we set 
the temperature T = 30 mK, a = 5 x 10~ 6 , and the unit 
of time is second (s). 

Fig. 2: The evolution of the off-diagonal elements of 
the reduced density matrix for the JCQ in the Ohmic 
bath. Here, we set B x = 51.8 fieV, B z = 0, T = 30 
mK, ujc = 5 (ps) _1 Hz, a — 5 x 10~ 6 , and the unit of 
time is picosecond (ps). The initial state of the qubit and 
environment see the text. 
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